arXiv:1502.07384v3 [cs.SI] 16 Sep 2015 


Spreading Processes over Socio-Technical Networks 
with Phase-Type Transmissions 

Masaki Ogura and Victor M. Preciado 


Abstract —Most theoretical tools available for the analysis 
of spreading processes over networks assume exponentially 
distributed transmission and recovery times. In practice, the 
empirical distribution of transmission times for many real 
spreading processes, such as the spread of web content through 
the Internet, are far from exponential. To bridge this gap be¬ 
tween theory and practice, we propose a methodology to model 
and analyze spreading processes with arbitrary transmission 
times using phase-type distributions. Phase-type distributions 
are a family of distributions that is dense in the set of positlve- 
valned distributions and can be nsed to approximate any given 
distributions. To illnstrate our methodology, we focus on a 
popular model of spreading over networks: the susceptible- 
infected-susceptible (SIS) networked model. In the standard 
version of this model, individuals informed about a piece 
of information transmit this piece to its neighbors at an 
exponential rate. In this paper, we extend this model to the case 
of transmission rates following a phase-type distribution. Using 
this extended model, we analyze the dynamics of the spread 
based on a vectorial representations of phase-type distributions. 
We illustrate our results by analyzing spreading processes over 
networks with transmission and recovery rates following a 
Welbull distribution. 


I. Introduction 

Understanding and controlling spreading processes over 
complex networks is an important problem with applica¬ 
tions in many relevant fields, including public health [1], 
malware spreading [2], and information propagation over 
socio-technical networks [3]. One fundamental result on the 
analysis of spreading processes over networks is the close 
connection between the spectral radius of the network and 
the dynamics of the spread [4]-[7]. Based on this result, the 
authors in [8]-[ll] proposed an optimization framework to 
find the optimal allocation of resources to control epidemic 
outbreak in different scenarios. 

The vast majority of spreading models over networks 
assume exponentially distributed transmission and recovery 
rates. In contrast, empirical observations indicate that most 
real-world spreading processes do not satisfy this assump¬ 
tion [12]-[14]. For example, the transmission rates of human 
immunodeficiency viruses present a distribution far from 
exponential [15]. In the context of socio-technical networks, 
the inter-arrival time of Twitter messages or the propagation 
time of news stories on a social media site can be explained 
using lognormal distributions [3], [12], [16], [17]. 
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There are only a few results available for analyzing 
spreading processes over networks with non-exponential 
transmission (and/or recovery) rates. The experimental study 
in [18] confirmed the drastic effect that non-exponential rates 
can have on the speed of spreading. In [19], an analytically 
solvable (although rather simplistic) model of spreading 
with non-exponential rates was proposed. An approximate 
analysis of spreading processes over graphs with general 
transmission and recovery times was proposed in [20] using 
asymptotic approximations. 

In this paper, we propose an alternative approach to 
analyze general transmission and recovery rates using 
phase-type distributions [21]. In contrast with [20], we pro¬ 
vide an analysis of general infection and recovery times for 
phase-type distributions without relying on asymptotic ap¬ 
proximations. The class of phase-type distributions is dense 
in the space of positive-valued distributions [22], hence, we 
can theoretically analyze arbitrary transmission and recovery 
rates. Furthermore, there are efficient algorithms to compute 
the parameters of a phase-type distributions to approximate 
any given distribution [21]. To validate our approach, we 
verify that the approximations in [20] are valid under certain 
irreducibility assumptions. The key tool used in our deriva¬ 
tions is a vectorial representations proposed in [23], which 
we use to represent phase-type distributions. 

The paper is organized as follows. In Section |I^ we state 
the spreading model under consideration. In Section [11^ we 
analyze this model when the recovery times follow a phase- 
type distribution, while the transmission times follow an ex¬ 
ponential distribution. Section IV is devoted to the dual case 
when the transmission times follow a phase-type distribution, 
while the recovery times follow an exponential distribution. 
Numerical simulations are presented in Section [V| 


A. Mathematical Preliminaries 

An undirected graph is a pair Q = (V,f), where V = 
n} is the set of nodes, and £ is the set of edges, 
consisting of distinct and unordered pairs {i,j} for i,j G V. 
We say that i is a neighbor of j (or that i and j are adjacent) 
if The adjacency matrix A G of ^ is defined 

as the jo, 1 {-matrix whose {i,j) entry is one if i and j are 
adjacent, 0 otherwise. Finally, the expectation of a random 
variable is denoted by £■[■]. 

We let I and O denote the identity and zero matrices with 
appropriate dimensions. Let e? denote the /th standard unit 
vector in and define Efj = ef{e’^)~^. By Ip we denote the 
p-vector whose entries are all one. We omit the dimension p 
when it is obvious from the context. A real matrix A, or a 



vector as its special case, is said to be nonnegative (positive), 
denoted by A > 0 (A > 0, respectively), if A is nonnegative 
(positive, respectively) entry-wise. The notations A < 0 and 
A < 0 are understood in the obvious manner. We denote the 
Kronecker product of matrices A and B by A ® B. Let A be 
a square matrix. The maximum real part of the eigenvalues 
of A is denoted by BiA). We say that A is Hurwitz stable if 
T 7 (A) < 0. Also, we say that A is Metzler if the off-diagonal 
entries of A are all non-negative. Finally, A is said to be 
irreducible if no similarity transformation by a permutation 
matrix makes A into block upper triangular matrix. 

Below we state some basic lemmas about Metzler matri¬ 
ces. The first lemma about the Hurwitz stability of Metzler 
matrices is standard and thus its proof is omitted. 

Lemma 1.1 ([24]): For a Metzler matrix A, the following 
conditions are equivalent. 

1) A is Hurwitz stable. 

2) There exists a positive vector v such that Av < 0. 

3) A is nonsingular and A^' < 0. 

The next lemma characterizes the marginal case when 
tj(A) =0. 

Lemma 1.2: Let A be an irreducible Metzler matrix. Then 
T 7 (A) = 0 if and only if A has a positive null vector. 

Proof: In the proof we denote the spectral radius of a 
matrix by p(-). Notice that, since A is Metzler, there exists 
c > 0 such that B — A + cI is positive and r\(B) ~ p{B). First 
assume that 7] (A) =0. Then the positive and irreducible ma¬ 
trix B satisfies p(B) = c. Therefore, by the Perron-Frobenius 
theorem (see, e.g., [25]), there exists a positive vector v such 
that Bv = cv, which implies Av = 0. Hence v is a positive null 
vector of A. The other direction can be proved in a similar 
way and hence it is omitted. ■ 

Finally we state the following lemma. 

Lemma 1.3: Let A,B G K”^". If A is irreducible and Met¬ 
zler and also B 7 ^ 0 is non-negative, then t 7 (A) < tj(A -|-B). 

Proof: The inequality holds true if A is non-negative 
by [25, Theorem 4.6]. If A is not non-negative, we consider 
instead A + cI where c > 0 is a scalar such that A + cI is 
nonnegative and also T7(A-|-c/) ~ p (A + cI). The details of 
the proof are omitted. ■ 

H. SIS Model with Phase-type Transmission 
AND Recovery Times 

In this paper, we investigate a popular model of spread¬ 
ing over networks called the susceptible-infected-susceptible 
(SIS) model. We provide a definition of this model as a 
family of continuous-time stochastic process below: 

Definition 2.1: Let Q be an undirected graph with n nodes. 
We say that a stochastic process (* = 1, • ■ ■ ,n) tak¬ 

ing values in the set {susceptible,infected} is a susceptible- 
infected-susceptible (SIS) model over Q if the process satis¬ 
fies the following conditions: 

Dl) For all f > 0 and i G there exists a ran¬ 

dom number R^'^t) > 0 such that, once z, becomes 
infected at time f, it remains infected during the 
time interval and becomes susceptible at 

time t-\-R^‘'> (t). We call /?(')(?) a recovery time. 


D2) For all t and /, there exists a renewal process 
0 = (t) < t']''' (?)<••• such that, for every neighbor 

j of i and k > 1 , if then zj becomes 

infected at time t-\-T^‘\t). We call the numbers T^‘\t) 
(k > 1) transmission times. 

D3) There exists a subset A C such that z,(0) 

is infected if / G A and z,(0) is susceptible otherwise. 
We regard that these initially infected nodes become 
infected at time 0, i.e., a node i G A has the recovery 
time R^‘\0) and the transmission times T^‘\q). 

We say that the infection-free equilibrium pft) = 
Pr(i is infected at time f) = 0 of the SIS model is exponen¬ 
tially stable if, as t ^ 0 °, pft) (i = 1,... ,n) converges to 
zero exponentially fastQ 

Throughout the paper we assume that all the recovery 
times and transmission times are independent to each other 
and, moreover, depend neither on t nor i. Therefore, we 
hereafter omit t and i and write R^‘\t) and T^‘\t) as R and 
Th when no confusion arises. 

In [20], the authors present the following two criteria 
for the stability of the infection-free equilibrium based on 
an asymptotic argument. First they show that, under the 
assumption that the inter-renewal times of T follow an 
exponential distribution with mean l/fi, if 

ij(A)<l/(i3£[B]), (1) 

then the meta-stable state of the infection probabilities must 
be equal to zero. They then show that, under the assumption 
that R follows an exponential distribution with mean 1 /5, if 

17 (A) < (1-/(5))//(5), (2) 

where / is the Laplace transform of the probability density 
function of the inter-renewal times of T, then the meta-stable 
state of the infection probabilities must be equal to zero. 
For the precise meaning of the meta-stability, the readers are 
referred to [ 20 ]. 

One of the primary purposes of this paper is to give 
justifications for the above claims, without asymptotic ap¬ 
proximations, under the assumption that either transmission 
or recovery times follow a phase-type distribution [ 21 ] intro¬ 
duced below. Consider a time-homogeneous Markov process 
in continuous-time with p + 1 (p > 1 ) states such that the 
states 1 , ..., p are transient and the state p -f 1 is absorbing. 
The infinitesimal generator of the process is then necessarily 
of the form 


where S is an invertible Metzler matrix with non¬ 
positive row-sums. Let [^] G {(j) G M.P) denote the 

initial distribution of the Markov process. Then, the time 
to absorption into the state p-fl, denoted by (0,5), is 
called a phase-type distribution. It is known that the set 
of phase-type distributions is dense in the set of positive 

*For simplicity in our presentation, we often say that i is infected (or 
susceptible) at time t when Zi(t) — infected (susceptible, respectively). 




valued distributions [22]. Moreover, there are efficient fitting 
algorithms to approximate a given arbitrary distribution by 
a phase-type distribution [21]. 

In order to analyze the conditions 0 and 0 using 
phase-type distributions, we shall state the following two 
assumptions; 

Assumption 2.2: R follows a phase-type distribution ((^,5) 
and the inter-renewal times of T follow an exponential 
distribution with mean l/j3. 

Assumption 2.3: The inter-renewal times of T follow a 
phase-type distribution {4>,S) and R follows an exponential 
distribution with mean 1/5. 


A. Vector Representations of Phase-type Distributions 

Vector representations of phase-type distributions play a 
crucial role in our analysis of SIS models. In order to 
introduce these representations, we first recall a vector repre¬ 
sentation of time-homogeneous Markov processes introduced 
in [23]; 

Lemma 2.4 ([23, Section II]): Let Q be the infinitesimal 
generator of a time-homogeneous Markov process taking its 
values in {l,...,p}. For distinct i and j let Nq^j denote 
the Poisson counter of rate qjj. Consider the stochastic 
differential equation 

p 

dx= Y, {Eml-Eu)xdNqi,^ 

i,m=l 


with the initial state a:(0) being a standard unit vector 
in K'’ with probability one. Then x is the time-homogeneous 
Markov process with the infinitesimal generator Q and the 
state space {ei,.. .,ep} C W. 

Using this lemma, we now provide a vector representation 
of phase-type distributions, as follows; 

Lemma 2.5: Let ((^,5) be a phase-type distribution. Let 
e^ denote the probability distribution on the standard unit 
vectors in M.P given by Pr(ei|i = e]) = (j)i for i= l,...,p. 
Consider the stochastic differential equation 

p p 

dx= Y (Emt - Eu)xdNsi,„, - Y EuxdNvf ( 4 ) 

i.m=l (.= 1 

with the initial condition that x(0) follows ei|i. Then, the 
extinction time random variable 


R = infjf > 0 ; x(f) = 0} 

= infjf > 0 ; 3t ,x{t) = e\ and jumps at time f} 


(5) 


follows ((p,S). 

Proof: The second identity in Q can be checked from 
the differential equation Q. Let us show the first identity. 
By Lemma 2.4 the stochastic differential equation 


dx= Y (Kt -E'ii 

i,m=\ 


p+1 


)x5Vs,„ -b Y (Ktle - 

i=\ 


represents a time-homogeneous Markov process with state 
space {ei,... ,ep+i} C and the infinitesimal generator 


in 0. Therefore, identifying the absorbing state, with 
the zero vector, we can see that the equation 

d^= YKT-E!^')xdNs,,,, + Yi-E’;t')^dK, ( 6 ) 

i,m=\ £=1 


represents a time-homogeneous Markov process with state 
space {ei,... ,ep,0} C and the infinitesimal generator 
in 0. Since in Q the last element of the variable x plays no 
role, the equation 0 is equivalent to the stochastic differ¬ 
ential equation Q with the state space Thus, Q gives 
the time-homogeneous Markov process with the infinitesimal 
generator 0 and the state space {ei,.. .,ep,Q} C M.P, where 
0 G is the absorbing state. Since x(0) follows e,i,, by the 
definition of phase-type distributions as the exit time, the 
random variable R follows (0,5). ■ 

Based on the above proved lemma, we can further give 
a vector representation of renewal processes whose inter¬ 
renewal times have a phase-type distribution. 

Lemma 2.6: Let (0,5) be a phase-type distribution. Let 
be the stochastic process that follows e^ at every time f > 0 
independently. Consider the stochastic differential equation 

p p 

dx^ Y {Emt- Eu)xdNsi,„ + Yi^<l>^J - Eu)xdNvg (7) 

i^m=\ t=\ 


with the initial condition that x(0) follows e^. Define To = 0 
and let 0 < Ti < 72 < • • • be the times at which x^ = 1 
and the counter Nvg jumps for some £. Then, the stochastic 
process T = {7]t}”^Q is the renewal process whose inter¬ 
renewal times follow (0,5). 

Proof: Let us first show T\ follows (0,5). By the defi¬ 
nition of 7j, on the interval [0,7’i), the stochastic differential 
equation (0 is equivalent to Q. In this equivalent equation, 
the random variables T\ and R are equal by their definitions 
and, furthermore, R follows (0,5) by Lemma 2.5 Therefore 
T\ follows (0,5). 

At time Ti, the differential equation Q reads dx = S(p—x. 
This means that x(7j) follows e^. Therefore, by using the 
memoryless property of Poisson counters, we can show that 
72 — 7j follows (0,5). In this way, an inductive argument 
proves that 7’ is a renewal process with its inter-renewal times 
following (0,5). ■ 


III. Exponential Transmission Times 


In this section we analyze SIS models under Assump¬ 
tion 2.2 and give sufficient conditions to achieve the expo¬ 


nential stability of the infection-free equilibrium. We notice 
that, under this assumption, D2 is equivalent to the following 
condition; 


D2') Whenever i and j are adjacent, i is susceptible, and 
j is infected, the node i becomes infected with an 
instantaneous rate of j3. 

The next proposition gives a vector representation of the 
SIS model under consideration. Throughout the paper, for 
each I <i<n and a real number A > 0, we let denote the 
Poisson counter with rate A. We assume that all the Poisson 
counters are independent of each other. 





Proposition 3.1: Let x(') (i = be the solutions of 

the stochastic differential equation; 


p 


Lm=l 


- Eux^‘'> dN^} + £0 (1 - f x(')) £ fly f dN^^'> 

t=\ j=\ 


( 8 ) 


with initial conditions: 

J x^'^(O) follows e^, i G A, 

I (0) = 0, otherwise. 

Define the stochastic processes Zi (/=!,...,«) by 


Zi{t) 


infected, IL''x*^')(f) = 1, 
susceptible, fx(')(f)=0. 


(9) 


( 10 ) 


Then, the processes z, are a SIS model satisfying Assump¬ 
tion 12,21 

Proof: Let fo > 0 and i be arbitrary. First assume that 
i is susceptible at time fo, i.e., fx(9(fo) =0. Then, from 
equation we see that x(9(f) = 0 at least until any of the 

counters jumps for some j such that ajj = 1. Until that 

time instant, the differential equation ([^ reads 

" ! -I 

£/x(') = £0 ^ Qij f xO') ^ (J1) 

f=l 


To this inequality we multiply from the left and obtain 
(f(l^x(9) = because Stp follows and 

is one of the standard unit vectors with probability one. 
Therefore, if a/j = 1 and fxl-'l = 1, then the quantity l^x(9 
becomes one whenever jumps. In other words, whenever 
i is susceptible, j is adjacent to /, and j is infected, the node i 
becomes infected with the constant rate of j3. Therefore, D2' 
is satisfied. 

Assume that i becomes infected at time fo. If fo = 0, 
then x(9(fo) follows by the initial condition (|^. On 
the other hand, if fo > 0, then the infection must occur 
by the transmission from a neighboring node and, by the 
argument in the last paragraph, such a transmission occurs 
when one of its neighbors j is infected and jumps. Since 
multiple Poisson counters cannot jump at the same time with 
probability one, the equation GD implies thatx(9(fo) follows 
Therefore, x^'^to) follows e(i) whatever value fo takes and 
hence we can, without loss of generality, assume that to = 0. 
Until i becomes susceptible, that is, until x(9 becomes 0, the 
stochastic differential equation ([^ is indeed equivalent to the 
vector representation 0 of the distribution ((j>,S), which has 
0 € as its exit state. This argument shows that the length 
of time until i becomes susceptible follows (0,5). Therefore 
D1 is also satisfied. 

Finally, the initial conditions 0 ensure that D3 is also 
satisfied. This completes the proof. ■ 

Using the stochastic differential equations 0 we can 
derive the following sufficient condition for the exponential 
stability of the infection-free equilibrium. 


Theorem 3.2: Under Assumption |2.2| if the matrix 

(0f) 

is Hurwitz stable, then the infection-free equilibrium of the 
SIS model is exponentially stable. 

Proof: Define = £’[x(9(f)]. Since pft) = 

it is sufficient to show that ^^‘\t) converges to 
zero exponentially fast for every i. Taking expectations in the 
differential equation ([^ (for details, see, e.g., [26]) yields 

— — = 52 “ Eu)Sim^ ^9 _ ^ Euvi^ 

t.m=\ £=1 


dt 


-f 

f=i 


(J)_ 


£a,7j3£[lV)lV^')]. 


Ignoring the last negative term in this equation, we obtain 
the inequality jdt < 5 "''^ (9 -|-j3 (A, ® (01'"))^, where the 
-valued function ^ is obtained by stacking ..., 

Therefore, we see that d^/dt < Ap^. Hence, if Ap 
is Hurwitz stable, then the comparison principle [27] shows 
that ^(f) converges to zero exponentially fast as t ^ as 
desired. ■ 

Using Theorem |3.2[ we can prove the validity of the 
condition Q under irreducibility conditions. 


Theorem 3.3: In addition to Assumption 2.2 assume that 


1.3 


the matrices A and S are irreducible and 0 is positive. Then, 
the condition Q is sufficient for the exponential stability of 
the infection-free equilibrium. 

Proof: First, we shall see that rj{Ap) is strictly in¬ 
creasing with respect to j3. In fact, for an arbitrary £ > 0, 
we have = Aj^ -|-£A(g) (01^). In this decomposition, 

Ap is Metzler and irreducible by the assumption. Moreover 
£A(g)(01^) is nonzero and nonnegative. Thus, Lemma 
shows ri{Ap)< ri{Ap+e)- 

Therefore, to prove the given claim, it is sufficient to 
show that ri(Apg) = 0 where ]3o = 1 /(£’[/?]tj(A)). This is 
equivalent to the existence of a positive null vector for Apg 
by Lemma [L^ because A is Metzler and irreducible. In the 
rest of the proof, we shall show that z = m 0 (—(5^^)^0) 
is such a null vector, where u is the eigenvector of A 
corresponding to the eigenvalue ri{A). Let us first show that 
z is positive. The vector u can be taken to be positive by the 
Perron-Frobenius theory because A is irreducible (see [25]). 
Moreover, since S is Metzler and Hurwitz stable. Lemma [LT] 
shows that 5^* does not have a positive entry. Also 5^* 
clearly does not have a zero row. Therefore, since 0 is 
positive, the vector — (5^')^0 is also positive. Hence z is 
indeed positive. Now, let us compute the product Ap^z: 

= -M® 0 - (/3otj(A)m) ® (0f (5^^)"^0) 

=-(M®0)(l-|-]3oT7(A)f’"(5^^)"^0). 

Since the mean of the distribution (0,5) equals — 0^5^*1 
([21]), we have l''(5^*)'''0 = 0'''5^'l = —Therefore, 
by the definition of j3o, the equation ( [T^ shows Ap^z = 0, 
as desired. This completes the proof of the theorem. ■ 


Remark 3.4: An advantage of Theorem 3.2 over the con¬ 
dition Q is that the theorem explicitly gives an upper bound 






on the speed of convergence to the infection-free state as the 
maximum real part of the eigenvalues of Ap. This would 
enable us to, for example, design the optimal strategies 
for distributing preventive resources over networks under 
constraints on the speed and the total amount of the resources 
available, as in [ 11 ]. 

IV. Exponential Recovery Times 
As the dual of the previous section, in this section we 
analyze SIS models under Assumption |2.3| Under this as¬ 
sumption, D 1 is equivalent to the following condition; 

Dl') Whenever a node becomes infected, it will recover with 
the instantaneous rate of 5. 

The next proposition gives a vector representation of the 


corresponding SIS model as in Proposition 3.1 

Proposition 4.1: Let (i = be the solutions of 

the stochastic differential equation 

ufx(') = -x(') dN^'’ + f; - Eu)x^^ dN^l 

^,m=l 

P n p 


£(e^ej -£a)x« - iV')) dN\ 


i=\ 


j=\ t=i 


(13) 

with the initial conditions 0 . Then the stochastic processes 
zu ..., Zn defined by are a SIS model satisfying 


Assumption 2.3 

Proof: Let to > 0 and i be arbitrary. Lirst as¬ 

sume that i is susceptible at time to. After time to, and 
while i is susceptible, the differential equation ( |T3] l reads 
= £01^=1 fl,; It ix[' 


dNij\ Multiplying l'" from the 


left yields d{t^x^‘^) = I^i■ Therefore, if 


Oij = 1 and xt = 1 , then, whenever the counter Nf/ jumps, 
the quantity becomes one, that is, i becomes infected. 

We then consider the case that i becomes infected at 
time tQ. As in the proof of Theorem 3.2 without loss of 


(;■) 


generality we can assume to — 0. After time 0, and while i 
is infected, the differential equation ( [T3] l reads 
p 


dx^^ =-x^^ dNf + £ {E,ni-Eu)x^^ dNll + 

i,m=\ 

£ {e^ej - Eie)x^‘^ dm’. 

e=\ 

By the first term of this equation, we see that x^'l becomes 
zero when and only when the counter jumps. This im¬ 
plies that i recovers with a rate of 5 and hence shows Dl' to 
be true. On the other hand, until N^g jumps, the variable x^'^ 
follows the same diffe renti al equation as Q. Since x(')(0) 


follows 60 , by Lemma 
Nv) jumps for some 


2.6 


the times at which x? = 1 and 


Lrom Proposition 
the exponential stability: 


Theorem 4.2: Under Assumption |2.3| if 

8>r]{l®S^ + {A + I)®{iiv^)), 


( 14 ) 


then the infection-free equilibrium of the SIS model is 
exponentially stable. 


Proof: As in the proof of Theorem 4.2 it is sufficient 
to show that ^^'\t) converges to zero exponentially fast as 
f oo for every i. Taking the expectation in ( [T3] l, we can see 
that 

d_^ 
dt 


l,m=\ t=\ 

+ ^ E ^ij - 01’" E E aijE[x^h^/\vi 

i=\ i=\ i=\l=\ 

Therefore d^/dt < (7® 5"'' -f (A 4-7) ® — dl)^. Hence, 

for pi{t) to converge to zero exponentially fast, it is sufficient 
that the matrix 7®5"'' 4- (A4-7) ®{<pv^) — 5l is Hurwitz stable 
by the same argument as in the proof of Theorem 3.2 This 
proves the sufficiency of the condition ( [l4| l. ■ 

Using Theore m |4.2[ we can then validate the effectiveness 
of the condition Qunder irreducibility conditions. 

Theorem 4.3: Under Assumption |2.3| if A is irreducible 
and V is positive, then (|^ implies the exponential stability 
of the infection-free equilibrium. 

Proof: Assume (|^. Then we have f{5) < 1 /(14-J 7 (A)). 
Since the probability density function of (0,5) has the form 
0^exp(5f)v (f > 0), we obtain 

())T(57-5)-^<l/(l4-t7(A)). (15) 

By Theorem |4.2[ it is sufficient to show that the Metzler 
matrix B = 7® 5 4- (A 4-7) ® (^v"'') — 5l is Hurwitz stable. 
Define z = u® ((57 —5)^*v), where m > 0 is the Perron- 
Lrobenius eigenvector of A. Since 5—57 is Hurwitz stable, 
the inverse (57—5)^* = —(5—57)^* is nonnegative and does 
not have a zero row by Lemma o Therefore the product 
(57 —5)^^v is positive because v is positive. Hence we see 
that z is positive. Now, using (HD), we can actually show that 
Bz < 0. Therefore, by Lemma B is Hurwitz stable. ■ 

V. Numerical Simulations 


(jj 

Vt 


t form the renewal process with its 
inter-renewal times following (0,5). This observation and 
the argument in the first paragraph of this proof prove that 
the stochastic processes z, satisfy D2. Also D3 holds true by 
the initial conditions (|9i. ■ 

we obtain the following criterion for 


In this section, we illustrate Theorem 4.2 through the 
comparison with the condition obtained in [20]. Let R follow 
a Weibull distributions with probability density function 
(a/5)(f/5)“^^exp(—(f/5)“), t > 0 , where a and b are 
positive parameters. In order to normalize the mean of the 
distributions to be one, we fix 5 = r(l-|-a^') where r(-) 
denotes the Gamma function. In [20], based on the condition 
it is concluded that the infimum of the recovery rate 5 such 
that the SIS model under Assumption |2.3| has the infection- 
free steady state equals Sq =r(l-|-a^^)r(a 4 -l)^/“t 7 (A)*/“. 

We compare the quantity 5o with the infimum re¬ 
covery rate 5i = T 7(7 ® 5"'' 4 - (A 4 - 7) ® (^v"'')), which 
we can obtain from Theorem |4.2| Let ^ be a real¬ 
ization of the Erdos-Renyi graph with 500 nodes. We 
vary the parameter a of the Weibull distribution as 












t 


Fig. 1. Weibull distributions fitted by phase-type distributions for a = 1.5, 
2.5, 3.5, and 4.5. Solid: phase-type distributions. Dashed: Weibull distribu¬ 
tions. 



Fig. 2. Minimum recovery rates. Circle: 5i. Square: Sq. 


a — 0.5, 1, 4.5, and 5. The Weibull distributions 

are fitted with phase-type distributions using the expecta¬ 
tion-maximization algorithm proposed in [21] (and available 
at http: / / home . math. au. dk/asmus/pspapers . html). 
Some of the fitting results are shown in Fig. [3 In Fig. g 
we compare the two recovery rates 5o and 5i. The graph 
shows the coincidence of the two rates except at a = 0.5. 

VI. Conclusion 

In this paper we have analyzed SIS models of spreading 
over networks with phase-type transmission and recovery 
times. We have derived sufficient conditions to tame the 
spread in terms of the eigenvalues of matrices that depend on 
both the graph structure and the parameters of the phase-type 
distribution. Our results mathematically justify the conditions 
found in [20] without using asymptotic arguments. The 
generality of the approach herein introduced is supported by 
the fact that the set of phase-type distributions is dense in the 
set of all the positive random variables. As a future work, we 
will develop control strategies to contain epidemic spreading 
with phase-type rates, as well as applications in the context 
of distribution of online content in social networks. 
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